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Abstract 

Multiple scattering and attenuation corrections in Deep Inelastic Neutron Scattering experiments are an- 
alyzed. The theoretical basis is stated, and a Monte Carlo procedure to perform the calculation is presented. 
The results are compared with experimental data. The importance of the accuracy in the description of the 
experimental parameters is tested, and the implications of the present results on the data analysis procedures 
is examined. 
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I. INTRODUCTION 



Since its creation in 1966 by Hohenberg and Platzmann , the Deep Inelastic Neutron Scattering 
(DINS) technique was considered the most direct probe of the momentum distribution of nuclei in 
condensed matter. The interest on this technique was further stimulated by subsequent develop- 
ments, that showed that features attributed to the interference between the neutron and the proton 
could be observed in experiments made on hydrogen 2 , thus determining the wave function of the 
protons directly from the experimental data. The availability of this technique as a customary re- 
search tool, opened in the last decade a new field for the investigation in condensed matter, and 
stirred up particular interest in the study of hydrogen dynamics, a topic for which this technique is 
particularly suitable. Despite the main activity in the field is held at Rutherford Appleton Labora- 
tory (United Kingdom), contributions from different laboratories were also reported in the pasl 3 - 4 - 5 - 
and recently 6,7,8,9 . 



The technique basically consists in an energy analysis carried out through the use of neutron 
resonant filters in the range of a few electron- Volts, and it is based on pulsed neutron sources. 
The spectra are recorded in time channels, which for the purpose of the study of momentum 
distributions must be translated to a momentum variable, thus obtaining the so-called Neutron 
Compton Profile (NCP). 

The theoretical basis of the technique was established by Sears 10 , who outlined the general 
procedure that must be employed to arrive to the NCP from the experimental data. Later, Mayers^ 
established the conditions of validity of the usually employed approximations. 

Recently, we re-examined the usual procedure to analyze the experimental data obtained from 
this techniqu e) 12 ' 13 , and and suggested improvements in the analysis of experiments involving light 
nuclei 14 . In the cited reference, we showed that the method to obtain the momentum distributions 
is, in general, a non-trivial task that involves a thorough knowledge of the different components of 
the experimental setup. In particular, it is important to know the energy spectrum of the incident 
neutrons, the detector efficiency, as well as an accurate description of the filter total cross section. 
However, in the above mentioned work we did not examine the sample-dependent effects such as 
multiple scattering and attenuation, which must be previously accounted for, before any analysis 
is attempted. It is a very well-known fact that multiple scattering and attenuation effects can be 
important even if all the reasonable cautions are taken in the sample design, since some low-signal 
portion of the observed spectra could nevertheless be seriously affected 15 . 
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Multiple scattering corrections is a long-debated subject, and has been extensively treated in 
the literature. Following the pioneering works of Vineyard 16 , and Blech and Averbach 17 , Sears 
thoroughly established its theoretical basis 18 , stating the integrals that describe the contribution 
of the n-th scattering process to the observed spectrum. The complexity of its solution was suc- 
cessfully tackled by Copley 19 who devised a Monte Carlo code, suited to a particular experimental 
situation. However, Copley's scheme could not be easily adapted to different kinds of experiments. 
On the other hand, the common knowledge in neutron scattering technique teaches that specific 
procedures must be devised for each particular experimental situatio n 1 ^ 20 . 

Therefore, there is a primary need of a customary correction tool for multiple scattering and 
attenuation effects in DINS experiments, that was only recently fulfilled. The authors have already 
presented numerical results from a new Monte Carlo code compared with experimental results in 
multiple scattering corrections in DINS experiments (see Refs.^ 9 -), and recently Mayers et ah— 
introduced the details of a new Monte Carlo code, related with the experimental setup thoroughly 
described in Ref.— . The Monte Carlo procedure described in the mentioned work contains a 
series of assumptions that were carefully analyzed in several publications. In the first place it is 
considered that the neutron has a well-defined final energy corresponding to the maximum of the 
main absorption of the resonant filter. In a recent analysis 23 we showed that such distributions 
are far more complex than considering a single final energy, and depends on the time channel 
and the dynamics of the scattering species. The scattering angles in^l are generated considering a 
random isotropic distribution. The validity of such assumption was investigated in Ref. 24 showing 
that the results obtained from such approximation deviate significantly from the exact ones when 
considering incident neutrons in the epithermal range. Another approximation contained in the 
mentioned work the description of the total cross sections of the scattering system with a constant 
plus a '1/v' absorption scattering law, which can deviate significantly from the real behavior when 
considering molecular systems 25 . 

In this paper we present the fundamentals of our Monte Carlo procedure to account for mul- 
tiple scattering and attenuation effects. Account is taken of the experimental details such as the 
energy spectrum of the incident neutrons, the resonant filter transmission and the efficiency of 
the detectors bank. Regarding the sample, inelastic scattering is taken into account employing 
suitable models for each analyzed case. The proposed model can be easily introduced as a double- 
differential cross section, either analytical or numerically, so a complete description of the energy- 
transfers and scattering angles are obtained. Experimental results of samples of different sizes are 
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shown, and the present code is benchmarked. The importance of an accurate model to describe the 
neutron- sample interaction is stressed, and finally implications of multiple scattering corrections 
on different situations are discussed. 



H. BASIC FORMALISM 

A basic description of the kind of experiments that we will analyze can be found in Ref . so 
we will give only a brief account here. We will analyze the case a DINS experiment performed 
in a pulsed source with an inverse-geometry configuration, i.e. the resonant filter is placed in the 
path of the neutrons emerging from the sample. A typical experiment consists of alternative 'filter 
out' and 'filter in' measurements, whose difference gives the NCP in the time- of- flight scale. 

We define Eq as the incident neutron energy (characterized by a spectrum 4>(£'o)), E its final 
energy, Lq the source-sample distance, L s the sample-detector distance and the scattering angle. 
The total time elapsed since the neutron is emitted from the source until it is detected, for a single- 
scattering event, is 

where m is the neutron mass. 

The resonant filter will be characterized by a total cross section <3f(E), a number density n 

and a thickness T, so the fraction of neutrons transmitted by it will be exp (— nT<5p(E)). If 
d 2 o 

(Eq,E,Q) is the sample double-differential cross section, then the difference count rate 

('filter-out' minus 'filter-in') at time of flight t, will b o 12 ' 26 



d^a 

c(t)= / dE ®(Eo)——(E ,E,Q)e(E)(l-e-» T ^) 
JE in f dildE 



-inf 
t=const 



dE 



dt 



AQ, (2) 



where £(£) is the detector efficiency, and AQ. the solid angle subtended by the detector. 

Integral © must be calculated at constant time, taking into account relationship ([TJ between Eq 
and E, and the Jacobian | dE /dt | must be evaluated at a fixed time 26 . The lower limit of integration 
is determined by the condition that in the second flight path the neutron has infinite velocity, i.e. 
E- m f = \mL\jt 1 . It is worth remarking that Eq. is a valid expression only if single scattering 
events would take place. 
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III. MULTIPLE SCATTERING 



In this section we will give an outline of the basic equations which govern the multiple scatter- 
ing processes of n-th order, that are employed in the Monte Carlo programs. For a more detailed 
treatment the reader is referred tc 15 and 18 . 

We will suppose throughout this paper an isotropic sample. Let S(<2, co) be the scattering law of 
the sample, Eq and E the incident and final neutron energies (being ko and k their corresponding 
wave vectors), dQ. the element of solid angle in the direction of the scattered neutron, the 
bound-atom scattering cross section of the sample (considered monatomic), and N the number of 
scattering centers, then the double-differential cross section is 

d 2 <5 N(5b k , 

defined as the average number of scattered neutrons with final energies between E and E + dE, and 
within a solid angle d£l, per unit incident flux. As usual we define Q = ko — k and /zco — Eq—E. 
The integral of Eq. © over all angles and final energies gives the microscopic total cross section 
a(£o). 

The above definition corresponds to the ideal textbook case where there is not multiple scat- 
tering. Turning to the real case let us define the macroscopic double-differential cross section as 
the probability that an incident neutron with a wave-vector ko will emerge from the sample with a 
wave- vector k 18 . In this definition we do not take into account neutrons non-interacting with the 
sample (i.e. transmitted). Its expression thus reads 

= s(k ,k), (4) 

dQ.dE 4nAk v °' ; ' 

where A is the cross-sectional area perpendicular to the incident beam. s(ko,k) is an effective 
scattering function that admits a decomposition in a part due to singly- scattered neutrons in the 
sample si(ko,k), another due to singly- scattered neutrons in the container sc(ko,k), and a third 
due to multiply scattered neutrons (with any combination of sample-container scattering events) 
SAf(ko,k) 

s(ko,k) =5i(ko,k)+5 M (ko,k)+5 C (ko,k). (5) 

The single scattering component si is simply related with the scattering law through the relation- 
ship 

Sl (ko,k)=N<5 b S(Q,e)H(ko,k), (6) 



where i/(ko, k) is the first-order attenuation factor, defined as the fraction of single-scattered neu- 
trons that fail to leave the sample due to multiple scattering and nuclear absorption 18 or that are 
not detected due to the detector efficiency 15 . Expression © inserted into ©, gives the real NCP 
including multiple scattering components. Its calculation, will normally involve a numerical sim- 
ulation based on the Monte Carlo method. 



IV. MONTE CARLO CODE 



In this section we will describe the numerical simulation devised for DINS experiments. Its 
fundamentals are based on Copley's method 19 , and they are extensively developed in 27 and 28 . 



A. Neutron Histories 



Neutron histories are generated with an initial unity weight. The incident neutron energy is 
decided randomly using the experimental neutron spectrum as the probability distribution. The 
flight path x for a neutron with energy E is given by the probability 

L t (E,x)T(E,x) 

P[h ' X) l-*(E,d) ' U) 
where the probability has been biased so the neutron never gets out of the sample 27 . In Eq. ©, 
T, t (E,x) is the macroscopic total cross section of the sample-container set a distance x away from 
the neutron previous scattering position, taken in the current flight direction, T (E,x) is the fraction 
of noninteracting (transmitted) neutrons in that direction after traversing a distance x, and d is the 
distance to the sample surface in that direction. To compensate the bias in the probability, a weight 
is assigned to each neutron which decreases according to the transmitted fraction in the traversed 
path, being 1 the initial value. Given the weight at step i — 1 the weight at step i is calculated as 28 

W( - = W/ -i(l-^(2M))|7§4, (8) 

where T, s (E,0) and Z t (E,0) are the macroscopic scattering and total cross sections, respectively, 
at position i — 1 and its ratio indicates the probability that the neutron will not be absorbed in the 
considered path. A history is finished when the weight drops under a predetermined cut-off value, 
so the number of scattering events is not predetermined. 

The assignment of new energies and flight directions is made via the use of model distributions 
for the double-differential cross sections of the sample and the container environments, normalized 



with the total cross section at the current energy Er± 



Ei,E,Q) 



No b k 
4no(Ei) h 



S model (Q, ®)- 



(9) 



B. Scoring 



At each step, the contribution of the current history to the detectors is calculated for each time- 
of-flight channel t. The final energy E to be considered for this channel is obtained from 15 



where N is the number of scattering steps, and L, is the flight path of step i, which was covered 
with an energy Ej. 

The quantity to be scored is the current weight, times the transmission factor from the cur- 
rent position to the sample surface in the direction towards the detector position, times the filter 
absorption ratio, times the detector efficiency 



It can be shown that the average of n after a large number of histories is the sought solution of 
Eq.© for the case of a macroscopic sample (Eq.©). 

C. Summary of input data 

The above described Monte Carlo procedure, requires a detailed description of the experimental 
setup and total cross sections of the involved materials. It also makes use of models for the scat- 
tering laws to describe the sample and container interaction with neutrons. Here we summarize 
the input data needed to perform it. 

• Incident spectrum as a function of energy. 

• Total cross section of the sample and the container materials as a function of the energy. 
These data must be tabulated in an energy range wide enough to cover not only the energies 
corresponding to the incident spectrum (epithermal energies), but also to consider the energy 
transfers after a number multiple scattering steps (typically thermal energies). 




(10) 



Zi = WiP(E h E,e)<l(E,d)(l-e 



T ° F ^)e(E). 



(11) 
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• Mean free path of the sample and the container as a function of energy. 

• Detector bank efficiency as a function of energy. 

• Input parameters for the chosen models to describe the scattering law of the sample and the 
container. Alternatively the models can be defined through a numerical input. 

• Geometry parameters for the proposed experimental setup and sample environment. 

• Total cross section of the resonant filter in an energy range comprising thermal neutrons (to 
give a good description of the ' 1/v' region 36 ), to energies above the main resonance (to 
describe the lower time-of-flight channels). 

V. EXPERIMENTAL SETUP 

The experiments were performed at the Bariloche pulsed neutron source (Argentina). Neu- 
trons, produced by the interaction of the electrons accelerated by the LINAC on a lead target, are 
moderated in a 4 cm thick polyethylene slab. A cadmium sheet is placed in the incident beam, to 
absorb thermal neutrons. The LINAC was operated at a 100 Hz rate. A collimated neutron beam 
1 inch diameter was employed. 

A schematic view of the of the DINS facility is shown in Fig. [TJ A movable cylindrical 
indium filter 0.25 mm thick, is placed in the flight path of the scattered neutrons. The movement 
is controlled remotely to perform alternative 'filter-in' and 'filter-out' measurements every 10 
minutes. The detector bank consists of six 3 He proportional counters (10 atm filling pressure, 
6 inch active length, 1 inch diameter) placed at a mean scattering angle of 56°. The detectors were 
covered with cadmium cylinders to minimize the background due to thermal neutrons. The flight- 
path lengths were 504 cm (source-sample distance), and 27.5 cm (sample-detector) respectively. 

The incident spectrum, measured employing a 3 He detector placed perpendicularly to the direct 
beam, is shown in Fig. |2]where the detector efficiency effect was accounted for. The detector bank 
efficiency was determined through the ratio of the spectrum of scattered neutrons on a lead sample 

36 Here we refer to the slow neutron regime, where the absorption cross section is inversely proportional to the neutron 
velocity. See Refi^i. 
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FIG. 1: Experimental setup employed for DINS experiments. The detectors are placed at a scattering angle 
of 56°. 

(which is mostly an elastic scatterer), and the spectrum measured on the direct beam. The result 
is shown in Fig. |3] where the cutoff near 0.5 eV is due to the cadmium cylinders that cover the 
detectors. 



VI. RESULTS AND DISCUSSION 



In this section we will analyze our experimental and numerical results on the Compton profiles 
for coin-shaped graphite and polyethylene samples of different sizes at room temperature. In 
the further paragraphs we will show results on samples whose sizes were chosen to serve as a 
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FIG. 2: Incident neutron spectrum multiplied by the energy in order to represent it in logarithmic Eq scale. 
Dips due to resonances in the cadmium sheet are observed. 
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FIG. 3: Detector bank efficiency. The cutoff about 0.5 eV is due to the cadmium cylinders which cover the 
detectors. 
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FIG. 4: NCP for the two measured graphite samples. Normal line: single scattering component; dotted line: 
multiple scattering; thick line: total scattering. Insets: attenuation factors 

benchmark on the numerical simulations, and they are not intended to represent optimized choices 
in the experimental design. Finally we will show the importance of the present corrections in 
thin samples, suitable for the experimentalists' work. Numerical simulations were performed 
employing the above mentioned Monte Carlo code, making use of the experimental parameters 
mentioned in the previous section, and using total cross section data for the Indium filter from 
Ref. 30 . 

In Fig. H|we show the results for the graphite samples, 3.54 cm diameter and 1 and 2 cm thick- 
ness respectively. A typical error bar is indicated for each experimental dataset. The approximate 
difference in count-rate ('filter out' minus 'filter in') at the peak maximum was 4 counts every 
10000 LINAC pulses for the thin sample and 6 counts every 10000 LIN AC pulses for the thick 
one. Measurements were carried out in 4 million LINAC pulses for the thin and 2 million LINAC 
pulses for the thick sample. Numerical simulations were carried out using a gas model for the 
graphite. This is a good approach at epithermal energies like in this case, with the condition that 
the temperature must be replaced by an effective temperature that takes into account the phonon 
dynamics 32 . The resulting effective temperature is 61.2 meV calculated on the basis of a Debye 
temperature of 1860 K 3 ^. In Fig. [4] we show the Monte Carlo results for the single and multiple 
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FIG. 5: NCP for the two measured polyethylene samples. The same notation as in Fig. |4] applies. Insets: 
attenuation factors 



scattering components as well as the total one. Besides the main peak at about 315 /jsec two extra 
peaks at 200 and 130 /usee are observed, due to the resonances of 3.85 and 9.07 eV of Indium, 
respectively. A good agreement is observed between the calculation and the experimental data, 
showing that the multiple scattering component has a peak shape that is broader than the main 
peak and is centered at lower times of flight, thus contributing to a significant distortion in the 
observed total scattering. In the insets of Fig. |4]we show the attenuation factors (Eq. ©), that 
have to be applied to the observed profile, once multiple scattering is subtracted. 

In Fig. |5]we show the same results for the two samples of polyethylene 3 cm diameter and 1 and 
2 mm thickness respectively, where typical error bars are shown. The approximate count rate at 
the peak position was 5 counts every 10000 LINAC pulses for the thin sample and 7.5 counts every 
10000 LINAC pulses for the thick one. Measurements were carried out in 5 million LINAC pulses 
for the thin and 4 million LINAC pulses for the thick sample. For the Monte Carlo simulations 
we employed the Synthetic Model 25 with the parameters from Ref. 33 , which was successfully 
employed to describe different integral magnitudes of the double-differential cross section. The 
model adequately describes the interaction between the neutron and the sample in different energy 
regimes, tending naturally to the commonly employed impulse approximation in the epithermal 
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FIG. 6: Integrated intensity for single (dashed line), multiple (dotted line) and total scattering (full line) 
for different thicknesses of graphite and polyethylene samples, compared with the experimental data (black 
circles). 

region. Although the shape of the observed main peak due to hydrogen is less affected by multiple 
scattering effects than in the case of graphite, it must be noted the distorting effect due to the 
attenuation factor that varies a 25 % from 200 to 350 /jsec thus affecting significantly the long- 
times tail of the Compton profiles. 

The general trend of the multiple scattering behavior can be analyzed by observing the total in- 
tensities observed in the main peak of the Compton profiles as a function of the sample thickness. 
In Fig. |6]we show the integral intensity of the main peak for single, multiple and total scattering 
obtained from Monte Carlo simulations at several sample thicknesses of graphite and polyethy- 
lene. In the same graph we show the results obtained from our experimental data. In the case of 
polyethylene, we measured a third sample 4 mm thick, that is not included in Fig. |5J but is shown 
in Fig. |6j It is worth to mention that the results from the simulations were multiplied by a constant 
(the same value in all the cases) in order to fit the experimental data. In both systems, we observe 
that the trend of the peak intensity as a function of the sample thickness can be correctly accounted 
for only if multiple scattering processes are considered. 

It is worth to emphasize the importance of a good description of the filter total cross section and 
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FIG. 7: Main frame: Comparison of the Monte Carlo run for a 2mm thick polyethylene sample (thick line) 
with a black detector model (thin line) and with a Lorentzian shape for the filter total cross section (dotted 
line) and the experimental data (dots). Lower frame: detail of the multiple scattering components. Inset: 
Lorentzian approach (dotted line) compared with the real total cross section for the indium filter (full line). 

the detector efficiency. For that purpose we performed simulations assuming two different cases: 

(a) a black detector (i.e. a detector with unit efficiency) and the filter described by the complete 
absorption cross section; 

(b) a filter described with a Lorentzian shape 36 , with the real detector efficiency. 

In both cases commonly employed approaches are used. The results are shown in Fig . |7]and 
in the inset the Lorentzian used to represent the filter together with the exact cross section 3 ^. We 
observe that both assumptions are inadequate and affect both the single and the multiple scattering 
components. In the case (a) the defect is manifested in an inaccurate description of the long-times 
tail. On the other hand, in (b) we observe an incorrect description in the short time region and the 
long-time tail, is mostly unaffected. It is worth to discuss the reason for both behaviors. The long- 
times region of the observed profile is mainly composed by emerging slow neutrons, which are 
absorbed by the filter according to the above referred ' 1/v' behavior. In the case of a black detector 
(a) our system is sensitive to those neutrons, while in (b) and our detectors are covered with 
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cadmium cylinders, whereby our detection system is insensitive to such neutrons. Both behaviors 
are observed in the calculated curves in Fig.|7J while experimental data only marginally illustrate 
the effect due to experimental errors. 

Finally, it is worth to mention that in common practice, the experimentalist will choose an ad- 
equate sample size in order to minimize multiple scattering effects, while keeping an acceptable 
signal-to-noise ratio. To illustrate the multiple scattering and attenuation effects on samples suit- 
able to the experimentalist, we present in Fig. [8] the results of our Monte Carlo program for thin 
sample thicknesses of graphite and polyethylene of the same diameters as presented in Figs. |H and 
|5J In the upper frame we show our results for graphite (1 mm thick). We observe that although 
the multiple scattering contribution is small, it is non negligible, and it will have to be properly 
computed if accurate values of the peak-shape parameters are to be obtained from the experiment. 
Special attention must be payed to the attenuation factor, that still has an appreciable variation over 
the range of times of interest. The reason is that i/(ko,k) in Eq. © contains not only the attenua- 
tion in the sample (negligible for a thin sample) but also the fraction of detected neutrons (detector 
efficiency effects). In the case of polyethylene, we show the results of a 0.15 mm thick sample. 
The multiple scattering effect is barely visible, but the attenuation factor has also an appreciable 
variation over the range of interest, thus affecting the observed peak shape. 



VII. DISCUSSION AND CONCLUSIONS 



Throughout this paper we examined different aspects that affect multiple scattering and at- 
tenuation effects in DINS experiments. We presented a Monte Carlo procedure that adequately 
describes Multiple Scattering and attenuation processes in DINS experiments. To attain a good 
agreement between the numerical simulations and the experimental data, accurate descriptions of 
the incident neutron spectrum, the detector efficiency as a function of the energy and the filter total 
cross section were necessary. These considerations add up to those stated in Refs . 12it ? , regarding 
the need of a good description of the experimental setup and the inadequacy of the convolution 
approximation in the description of the neutron Compton profiles. It must be remarked the impor- 
tance of an accurate description of the efficiency of the detector system. As such, we understand 
the detecting setup that comprises the kind of detectors employed, the geometry and the materi- 
als involved. For instance, the detector system employed in this work is composed by 3 He tubes 
covered with cadmium, that result in the efficiency shown in Fig. The importance of the knowl- 
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FIG. 8: NCP for thin graphite and polyethylene samples chosen to approach the experimentalists' needs. 
The same notation as in Fig. @]applies. Insets: attenuation factors 

edge of the detector efficiency as a function of the energy can be understood in the light of the 
recent work on the analysis of final energy distributions in DINS experiments 2 As a result of this 
analysis it was concluded that at every time of flight a distribution of final energies is operative 
in the neutron Compton profiles, instead of a single well-defined energy as is usually assumed 
in the customary data processing procedures 21 . As a consequence is is not sufficient to assume a 
detector operative at a single energy (in which case the knowledge of the efficiency function would 
be irrelevant), but the behavior of the efficiency function is essential, as confirmed by the results 
presented in Fig. |7] 

The results presented on Multiple Scattering effects should be considered when designing a 
DINS experiment. In the analyzed graphite samples, the multiple scattering components have a 
significant structure that affect the shapes of the neutron Compton profiles. In the case of polyethy- 
lene, multiple scattering contributes with a flat background, but on the other hand the attenuation 
factor has a significant variation over the time range of interest. This is a consequence of the sub- 
stantial variation of the efficiency function of our detection system in the range of energies shown 
in Fig. |3] 
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From the analyzed examples, we conclude that a wise choice of the sample thickness is still a 
valid rule, taking into account that a good contrast between 'filter-out' and 'filter-in' positions is 
required. Even if the sample sizes are adequately chosen, multiple scattering and/or attenuation 
corrections will necessarily have to be considered. The case of Fig.[5]is illustrative. Even if the 
multiple scattering contribution is small, the correction due to the attenuation factor can still be 
important (as shown in our case), given that it also includes the detectors' efficiency correction^, 
which in our case has an appreciable variation with the energy. In general this result will depend 
on each particular detection system. All the considered cases show the need to perform accurate 
multiple scattering, attenuation and efficiency corrections. To this end numerical simulations are 
the most adequate procedure. In this paper we presented a suitable correction tool, for which an ex- 
perimental benchmark with samples considerably affected by these corrections was satisfactorily 
performed. 

These corrections will normally have to be taken into account before proceeding to the data 
analysis, i.e. obtaining the kinetic energy distributions of the atoms, peak areas, etc. Particular 
importance will have the corrections in the analysis of lighter nuclei. Monte Carlo simulations on 
heavier nuclei not shown in this paper 21 reveal that the multiple scattering component tends to be 
located below the main peak and it is roughly proportional to it, thus having negligible distortion 
effects on its shape. This consideration must be brought together with those mentioned in Ref. 12 
regarding the extreme care that must be exercised when analyzing light nuclei with the DINS 
technique. 
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